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Considerations on Neuberger's operator. 
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We discuss new approaches to the numerical implementation of Neuberger's operator for lattice fermions and 
the possible use of block spin transformations. 



Recently major progress has been made in the 
lattice discretization of fermionic fields. Starting 
from the overlap formulation, Neuberger has pro- 
posed a lattice Dirac operator which avoids the 
additive renormalization of fermion masses [y. 
This operator satisfies an identity, the Ginsparg- 
Wilson relation |g|, by which one can define ex- 
act lattice chiral symmetry Q . Other lattice op- 
erators which satisfy the Ginsparg- Wilson rela- 
tions have been obtained as implementations of 
the "perfect action" M . In this note we will con- 
centrate on Neuberger's operator and we will fol- 
low his notation. If one defines the Wilson lattice 
operator as 

V'(a:) 



(Dwi'Kx) 



-Y,[{^-%.)U,ixmx + ^i) 



+ (1 + l,,)Uf,{x 

Neuberger's operator is 
1 



fi)ilj{x~ n)] 



D 



-Sl + V) 



(1) 



(2) 



where the unitary operator V is given by 

V = Dw{dI,Dw)-^^^ (3) 

Provided the hopping parameter k in Eq. ^ 
is chosen within a suitable range (1/8 — 1/4), 
Neuberger's operator defines a single flavor of 
massless lattice fermions, symmetric under chiral 
transformations. 

Since in actual calculations the size of the ma- 
trix Dyy is very large, for a practical implemen- 
tation of Neuberger's operator it is crucial to find 
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computationally efhcient ways to calculate its ac- 
tion as well as the action of its inverse on any 
given vector. In some sense, the domain wall for- 
mulation of lattice fermions ^-0] , which preceded 
the introduction of the operator of Eq. |[ pro- 
vide a numerical procedure for its implementa- 
tion. Other approximate methods have been pro- 
posed in [p|-pl| . Most of these numerical methods 
proceed through an approximation to the opera- 
tor 



eiH) ^ H{H^) 



-1/2 



^^Dw{d\^D 



W^W) 



-1/2 



(4) 



where H = 75-Dh/ is a Hermitian operator. Here 
we would like to focus directly on the operator 
V . We will make the assumption that D\^Dv/ 
has no zero eigenvalues, so that {D\^Dw)~^^^ is 
well defined (it can be argued that such zeroes 
are exceptional and irrelevant for the continuum 
limit). 

Given any square matrix M such that IvP M 
has no zero eigenvalues, its "polar decomposition" 
is defined by 



M ^ UA 



(5) 



where the matrix U is unitary and the matrix A is 
positive-definite Hermitian |12| . It mirrors the ex- 
pression of a complex number in terms of its phase 
and modulus. The polar decomposition is unique, 
with A = {M'^Mf/^ and U = M{M'^M)-^''^ 
From Eq. ^ we see that the matrix V is the uni- 
tary factor in the polar decomposition of the Wil- 
son lattice operator Dw- The main observation 
which we would like to make is that the unitary 
factor U in the polar decomposition of a matrix 



M satisfies a maximum principle, namely U is 
the unique matrix which maximizes the expres- 
sion Re Tt{U'M^), where U' ranges over the en- 
tire space of unitary matrices of the same size as 
M: 



ReTi (UM^) = Max 



{U'U'^=I} 



Re Tr([/'Mt) (6) 



The proof is straightforward. We write M in its 
polar form and write U' as U' = UW, where W 
is also unitary. This gives 

Re Tt{U'M^) = Re TilUWiA'PAiy^^U^] 

= ReTr[W^(MtM)i/2] (7) 

But (M^MY^^ is a positive-definite Hermitian 
matrix and we can therefore choose a basis where 
(M'lM)^^^ = Diag(Ai,A2,...Ajv), A, > 0. In 
this basis Re Tr[H/(MtM)i/2] = J2^Re{w^^,)X, 
and since the eigenvalues Xi are all positive and 
the unitarity of W implies \wi^i\ < 1, the maxi- 
mum occurs for W — I, i.e. U' — U. Although 
this proof is very simple, the property expressed 
by Eq. M can hardly be found in books on ma- 
trix algebra jl3|. Indeed very few pages, if any, 
are generally devoted to the polar decomposition 
itself. 

Maximization and minimization principles 
have often proven very useful in numerical analy- 
sis, since they can form the starting point for ef- 
ficient schemes of approximation. Here we would 
like to briefly mention a few possibilities. Assum- 
ing that one is interested in calculating either Dx 
or D~^Xi one could construct the Krylov space 
spanned by the vectors obtained acting with H 
on X (this space also underlies the approximation 
methods of §-0): 



Xk = H'^x 



0. 



(8) 



In order to respect the symmetry properties of V 
under 75 transformation, it is desirable to aug- 
ment the basis by considering separately the two 
projections 



1 



Xfcs 



S75 



-Xk 



±1 



(9) 



For sufficiently large n this basis would become 
overcomplete, but we must work under the as- 
sumption that values of n much smaller than the 



dimensionality N of the full vector space will pro- 
duce reasonable approximations. We will finally 
denote by rjks an orthonormalized basis in the 
space spanned by the vectors Xks- (This will re- 
quire forming a number of scalar products of or- 
der n^ . We are working under the hypothesis that 
the values of n one must consider make it feasi- 
ble both to calculate these scalar products and 
to store the vectors r/^g.) At this point the max- 
imization of Re Tr([/'M'I') may be restricted to 
unitary operators in the space E spanned by the 
vectors "qks, namely 



U' = 



y^ Vks,k's''nksr]i' 
ks,k' s' 



(10) 



where /^ denotes the identity in the complement 
of E. Maximizing Re Tr(t/'Mt) one finds that 
the (2n) x (2n) matrix Vks,k's' is the unitary fac- 
tor in the polar decomposition of the truncated 
matrix r]l,, ^, DwVks ■ As already mentioned above, 
this procedure will work only if reasonably small 
values of n can produce satisfactory approxima- 
tions. 

The truncation of Eq. |o| does not represent, 
however, the only way to take advantage of the 
maximization principle. Another possibility is to 
assume an approximately ultralocal form for V. 
V is known to be local |10| and the approximation 
would consist in making it ultralocal. Imposing 
such a condition on a unitary operator may be 
problematic, but one could write V = expiiK) 
with K'^ — K and truncate to the matrix el- 
ements of K which exceed a maximum separa- 
tion. Otherwise, one could try to construct a 
better approximation to V by the refinement of 
a first approximation: V = V' V'^^\ A particu- 
larly appealing possibility is to take for y(°) the 
operator defined by a coarsening of the lattice 
obtained by a blocking similar to those used in 
multigrid algorithms. Multigrid techniques have 
been tried for lattice fermions, but their appli- 
cation to the Wilson lattice operator [Q turned 
out to be only marginally effective, most likely 
because of the presence of an additive mass renor- 
malization. Neuberger's operator is not affected 
by a similar renormalization and this offers the 
hope that multigrid methods may work much bet- 
ter. In a multigrid approximation, one divides the 



lattice in cells of 2'^ sites (for one level of coars- 
ening). The gauge must be fixed within the cells, 
in order to make the transport factors as close as 
possible to the identity. Since this gauge fixing is 
local to the cells, it is not computationally expen- 
sive. The local gauge fixing leaves the freedom of 
performing gauge transformations common to all 
the sites within the cells. This becomes a gauge 
symmetry of the blocked lattice. After the local 
gauge fixing it makes sense to define a projection 
over the average cell fermionic fields. We define 
new basis vectors 



(0) 



-d/2 



E 



mx+y 



(11) 



where the vectors 'qx form a local basis for the 
fermions (spin and color indices are left implicit), 
X and X denote the integer valued coordinates of 
the original and coarse lattices, respectively, and 
y are d-dimensional vectors with components that 
take value or 1. We can now define 



y'''-T.-'xWxW^''+i^ 



B(0) 



(12) 



XX' 



where Ie{o) denotes the identity in the comple- 
ment of the space spanned by ry-^ . The ma- 
trix elements v-^-^, can then be determined by 
maximizing ReTr(F'^°)D|y). Alternatively, one 
could maximize ReTr(y('')£)Oo)t)^ .^^j^ere d''^ de- 
notes a coarse lattice Wilson, obtained by defin- 
ing transport factors U)^ between neighbor- 
ing cells as the SU(i) matrices that maximize 
ReTr(C/jf°^''S'^t)^ ^^^^^ g^ denotes the sum of 
the transport factors U^ connecting the bound- 
ary sites of the cells in the original lattice. In 
both cases, the size of the matrix to determine is 
reduced by 2^* x 2^^. 

Of course the crucial question is whether these 
or similar computational procedures can provide 
a sufficiently accurate approximation to Neu- 
berger's operator without excessive usage of CP 
time and memory. We are currently working to 
test the approximations outlined above and hope 
to report on our results in the near future. At 
the same time, we hope that this communication 
may open new avenues in the quest for efficient 
algorithms for lattice fermions. 
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